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Our goal is to build an Electrocardiogram (ECG) that not only calculates 
the heart rate automatically, but can also detect other heart abnormalities as 
well. This requires more advanced analysis of the ECG Signal. There are 
several steps that need to be accomplished in order to achieve this goal, as 
outlined in the flowchart below. 

Overall Data Flowchart for an Automatic ECG 


Breaks down the different steps which need to be accomplished in 
order to accomplish ECG signal analysis 


Data acquisition and signal conditioning are covered in [missing_resource: 
m18956][missing_ resource: m18957] 


Before we go into more detail about how to build an ECG, it is helpful to 
understand how the ECG works and how to interpret the data you receive. 


Physiological Background of the ECG 


An EGG is a non-invasive diagnostic device to monitor the condition of the 
heart through its electrical activity. This signal is acquired through 
externally located electrodes that adhere to the skin. A simple, clinical lead 


placement uses three leads: left arm, right arm and left leg (Figure 2). The 
electrical activity versus time forms an electrograph and can be used to 
determine and diagnose heart abnormalities and arrhythmias. This principle 
is based on Einthoven’s law. 

Electrode Placement for Three-Lead ECG 


Lead 1 is 
attached to the 
right wrist, Lead 
2 to the left wrist 
and Lead 3 is 
attached to the 
right ankle as a 
ground electrode. 


All recorded electrical activity of the electrocardiogram corresponds to the 
net electrical current in the heart over time, depolarizing parts of the heart 
in sequence. The electrical impulse is initiated by the sinoatrial (SA) node. 
This causes the atria to contract and is evident on the ECG as the P wave. 
Next, there is a delay caused by the conduction of the impulse to the 
atrioventricular (AV) node such that the physical contraction of the atria 
have time to complete before the contraction of the ventricles. The QRS 
complex on the ECG is due to the depolarization of the ventricles, and 
occurs when the ventricles contract. Finally, the T wave on the ECG is due 
to the repolarization of the ventricles (Pflanzer, 2004). Therefore, each 
heartbeat corresponds to a pulse on the ECG beginning with each P wave 
and the ending with each T wave. The heart rate can be determined by 


determining the time it takes to complete one beat and is typically reported 
in beats per minute. Figure 3 demonstrates the characteristic shape of the 
waveform in a healthy patient. 
Typical ECG Signal 

R 


An example of the 
typical shape and 
location of the 
various 
components of an 
ECG signal. 


Background Information on Monitored Heart Conditions 


To perform automatic detection of an ECG signal, there needs to be 
something that clearly delineates a certain abnormality from other signals. 
Therefore, the signal processing selected is to detect ventricular 
hypertrophy and old myocardial infarctions. Luckily, these types of 
abnormalities are both very useful for doctors to diagnose a patient and 
have distinguishable ECG features. 


Ventricular hypertrophy is the enlargement of either of the ventricles. 
Left ventricular hypertrophy is particularly common in athletes as well as 
an indicator of hypertension. It is also used in the Framingham risk equation 
to predict future cardiac problems the patient may face (ECG 
Abnormalities, 2006). One of its characteristic ECG patterns is the inverted 
T wave (Figure 4). 


An Example of an ECG Signal with an Inverted T-wave 
reerecaettveiaeett its 


A myocardial infarction (a.k.a. heart attack) is caused by the complete 
blockage of one of the coronary arteries. The coronary artery is what 
supplies the heart muscle with blood. A blockage prevents blood from 
reaching the surrounding muscular tissue resulting in necrosis. This damage 
is permanent so the resulting ECG characteristic will remain with the 
patient. Therefore, the doctor can easily tell if a patient has had a heart 
attack in the past. The ECG of an old myocardial infarction is characterized 
by a significant Q wave (Figure 5). This means that the Q peak is unusually 
deep, usually with amplitude of about one-third that of the R peak (Dubin, 
2000). 

An Example of an ECG Signal with a Significant Q-wave 


By building an ECG that can automatically detect these and potentially 
other heart abnormalities, it will be easier for doctors to monitor multiple 
patients. This is especially important in third world countries where there 
are often far too many patients in one clinic than one doctor or nurse can 
adequately care for. 


Collecting and Filtering Live ECG Signal 


Collecting ECG Signal: Hardware 


Using the NI Elvis breadboard and data acquisition system, the ECG signal 
is collected from the three leads. Lead 1 is connected to the right arm, Lead 
2 is connected to the left wrist and Lead 3 is connected to the right ankle 
(see [missing_resource: m18955] 

ECG Signal Collection Circuit 


Lead 1 
Output: to 
Lead2 o Filters and 
LabVIEW 


ECG signal collection circuit adapted from that of Radio 
Locman (www.rlocman.ru). The AD620 is used as a 
differential amplifier with a gain of ~7 in order to combine 
the electrical signals from each of the leads into one easily 
readable signal. 


In NI LabVIEW, the Data Acquisition Assistant (DAQ Assistant) is used to 
collect the signal after preliminary band pass filtering. The data is sampled 
at arate of 1 kHz. Our LabVIEW VI is available here. 


Signal Conditioning 


Signal filtering is necessary to help isolate the frequencies found in the 
ECG signal from the noise. With a three lead system, the majority of the 
noise comes from the electrical activity in the muscles on the arm, or 
electromyography (EMG) noise. EMG signals are present in a wide 
frequency band which overlaps with the ECG signal in the lower 
frequencies. Therefore, with this set-up, it is impossible to completely 
remove EMG noise from the ECG signal. Therefore, it is helpful for the 
patient to relax and remain still while the data is being collected. In 
addition, 60 Hz noise is present from power line interference which also 
must be removed. 


Analog Band Pass Filter Design 


The first stage of filtering is an analog filter built onto the NI Elvis 
breadboard. It is a bandpass filter with cut-off frequencies of 0.5 and 150 
Hz. This will help eliminate the high frequency noise from the muscles 
before the signal is greatly amplified. 

Analog Bandpass Filter Circuit Diagram 


Low Pass Filter High Pass Filter 
Cutoff Frequency: Cutoff Frequency: 


150 Hz 0.5 Hz 
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Bode Plot — Analog Band Pass Filter 


Magnitude Response 


Digital Filters Using LabVIEW 


Once the signal has been acquired by the DAQ Assistant into LabVIEW, it 
is processed by two additional filters and amplification of 100 times. The 
first filter is a band stop filter between 55 and 65 Hz to eliminate power line 
interference. A third order Butterworth (IIR) was used to implement this 
filter because it is low order and has a good frequency response for this 


signal. 


Magnitude of the Frequency Response — Digital Band Stop Filter 


Magnitude Response 
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The second is a tenth order Butterworth low pass filter. The cutoff 
frequency of this filter is 80 Hz to further eliminate EMG noise. The ECG 
signal is located between 0.5 Hz and about 70-80 Hz depending on the 
individual. 
Magnitude of the Frequency Response — Digital Low Pass Filter 
Magnitude Response 
= 


T 


O-| 


i i i 1 ! ' ! ! ! 1 
0 SO 100 150 200 250 300 350 400 450 500 
Frequency 


The resulting signal is displayed in real time in the graphical user interface, 
also designed in LabVIEW. 


Algorithms for ECG Signal Analysis 


ECG Signal Analysis: Abnormality Detection 


The flowchart below breaks down the tasks needed to accomplish signal 
analysis in greater detail. Both LabVIEW and Matlab were used to 
accomplish these tasks. 


Matlab Implementation of Filter Banks Analysis 


Many scenarios deal with signals which contain specific energy 
distributions in the frequency domain. For example, with regard to the 
ECG, a significant proportion of the energy from the QRS complex extends 
to a frequency of 40 Hz, and even more if the Q, R, and S waves have very 
sharp morphologies. The P and T waves, in general, have a significant 
proportion of their energy only up to 10 Hz. Thus, a strategic approach to 
detecting heartbeats is to analyze different sub-bands of the ECG, rather 
than just the output of one filter which maximizes SNR of the QRS. 


In this filter bank analysis technique we use 5 sub-bands, each one has 
bandwidth 5Hz. The ECG signal is processed by those 5 sub-band filters 
and downsampled. The downsampled signal is 
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A variety of features which are indicative of the QRS complex can be 
designed by combining sub-bands of interest from | = 0, 1, ... , M-1. For 
example a sum-of-absolute values feature P1 can be computed using sub- 
bands 1, 2, and 3. 
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P1 has a value which corresponds to the energy in the frequency band [5.6, 
22.5] Hz. Similarly, P2 and P3 can be computed using sub-bands {1, 2, 3, 
4}, and {2, 3, 4}, respectively, and these values are proportional to the 
energy in their respective sub-bands. Heuristic beat detection logic can be 
used to incorporate some of the above features which are indicative of the 
QRS complex. 


Figure 4 gives an overview of the sequential levels in the beat detection 
algorithm. The goal of the detection algorithm is to maximize the number of 
true positives (TP’s), while keeping the number of false negatives (FN’s) 
and false positives (FP’s) to a minimum. Since it is not possible to arrive at 
this goal using one simple detector, multiple detectors with complementary 
FN’s and FP’s performances are simultaneously operated and the results of 
each fused together to arrive at an overall decision. The advantage of this 
strategy is that multiple features which are indicative of the QRS complex 
can be used to detect beats. 

Beat Detection Algorithm Schematic 
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A one-channel beat detection block is described in Figure 5. 
Schematic of One-Channel Beat Detection Block 
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The detection strength Ds of an incoming feature (e.g., P1, P2, P3) is 
determined by comparing with the signal and noise levels (SL and NL, 
respectively): 
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If a feature’s value is less than NL then Ds is limited at 0, and if it is above 
SL then DS is limited to one. When a feature has a greater than a specified 
threshold (preset between zero and one) it is classified as a signal peak and 
the signal history is updated with the feature’s value. If the feature has a 
smaller than the threshold it is classified as a noise peak and the noise 
history is updated with the feature’s value. 


The detail operations of each level can be found in the paper by V. Afonso 
et al (ECG beat detection using filter banks, IEEE Trans. Biomed. Eng. 46: 
1999). Also, download and run ECGmain.m to test out filter banks. 


LabVIEW Programming for Signal Processing and Interfacing 


LabVIEW is used to not only process and acquire the signal, but to also 
develop an user friendly interface that automatically alerts the user to any 
abnormalities or arrhythmias detected. This particular ECG is developed to 
detect left ventricular hypertrophy and old myocardial infarction. It will 
also automatically display the heart rate and whether or not the patient’s 
heart rate is within a healthy range. 

The Graphical User Interface 
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This view of the the GUI shows a normal, healthy ECG signal. The 
indicator lights (6) will change color depending on the diagnosis. 


Additionally, the green trace turns red if any abnormality or 
arrhythmia is detected. 


The user interface consists of several useful features. First, there is a text 
input for the patient’s name to avoid confusion in the clinic. Second, the 
ECG chart is automatically updated with the name and live signal data. 
Third, the user can control the data acquisition periods. The heart rate and 
any abnormalities are displayed in the patient report. The warning lights to 
the left of the patient report correspond to the current condition of the 
patient’s ECG. If everything is normal, the green light is on and the signal is 
green. If an arrhythmia is detected, the yellow light turns on and a warning 
message will appear in the patient report that indicates the type. Similarly, 
the orange light turns on if an abnormality is detected and a warning 
message appears in the report. For either of these, the ECG trace turns red 
in the chart window. A red light indicates that no recognizable signal is 
being detected. This can mean that the patient has flat-lined, that the circuit 
board is not powered or that the leads are not connected properly. As a 
result, the warning message will ask the user to check the lead and power 
connections. Finally, the form can be reset between patients. You can 
download our VI here to test it out. 


